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ABSTRACT 

In the context of AGN unification scheme rapid variabihty properties play an 
important role in understanding any intrinsic differences between sources in differ- 
ent classes. In this respect any clue based on spectral properties will be very useful 
toward understanding the mechanisms responsible for the origin of rapid small scale 
optical variations, or microvariability. Here we have added spectra of 46 radio-quiet 
quasars (RQQSOs) and Seyfert 1 galaxies to those of our previous sample of 37 such 
objects, all of which had been previously searched for microvariability. We took new 
optical spectra of 33 objects and obtained 13 others from the literature. Their II/3 and 
Mg II emission lines were carefully fit to determine line widths (FWHM) as well as 
equivalent widths (EW) due to the broad emission line components. The line widths 
were used to estimate black hole masses and Eddington ratios, £. Both EW and FWHM 
are anticorrelated with £. Nearly all trends were in agreement with our previous work, 
although the tendency for sources exhibiting microvariability to be of lower luminosity 
was not confirmed. Most importantly, this whole sample of EW distributions provides 
no evidence for the hypothesis that a weak jet component in the radio quiet AGNs is 
responsible for their microvariability. 

Key words: galaxies: active - quasars: emission lines - quasars: general 



1 INTRODUCTION 

Rapid small scale optical variations, or intra-night mi- 
crovariability is a well known characteristic of Active Galac- 
tic Nuclei (AGNs) but the processes causing the bulk of 
these microvariations is still a matter of debate. The vari- 
ability mechanisms in the radio-loud objects are widely be- 
lieved to be connected to conditions in the relativistic jet. 
However, it is still unclear if in the radio-quiet objects the 
nature of intra-night variability is different, or whether the 
faint, variable jet is also a dominant component for them. 
Czerny et al. (2008) have tried to understand the microvari- 
ation mechanism within the framework of following scenar- 
ios: (i) fluctuations from an accretion disk (e.g., Mangalam 
& Wiita 1993); (ii) irradiation of an accretion disc by a 
variable X-ray flux (e.g., Rokaki, CoUin-SoufTrin & Magnan 
1993; Gaskell 2006); and (iii) the presence of modestly mis- 
aligned jets in radio-quiet quasars or a "blazar component" 
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(e.g., Gopal-Krishna et al. 2003). They concluded that the 
blazar component model is the most promising to give rise 
to intra-night optical variability. 

In this blazar component scenario, spectral properties 
of the sources can play a crucial role in constraining the 
models further. For instance, if blazar components are dom- 
inating the variability of RQQSOs, then, due to the increase 
in the continuum produced by the jets one expects smaller 
equivalent widths (EW) of prominent emission lines such as 
H/3 and Mg ll for sources with microvaraibility as compared 
to their average value in a sample including non-variable 
sources. Recently, in Chand, Wiita & Gupta (2010; here- 
after Paper I), we have worked toward this goal by ex- 
ploiting the optical spectra available from the Sloan Digi- 
tal Sky Survey, Data Release 7 (SDSS DR-7; Abazajian et 
al. 2009). We carried out careful spectral modeling of the 
H/3 and Mg ll emission line regions for RQQSO and Seyfert 
1 samples already searched for microvariability (hereafter 
referred to collectively as RQQSOs). In Paper I, we first in- 
vestigated any effect of key spectral parameters (e.g., EW 
and FWHM) on the microvariability of RQQSOs. Second, 
we estimated other relevant AGN parameters such as the 
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black hole mass, M;,^ and Eddington ratios. To do this, 
their H/3 and Mg ll emission lines were carefully fit to de- 
termine line widths (FWHM) as well as equivalent widths 
(EW) due to the broad emission line components. The line 
widths were used to estimate black hole masses and Edding- 
ton ratios. The EW distributions did not provide evidence 
for the hypothesis that a weak jet component in the RQQ- 
SOs is responsible for their microvariability, and perhaps 
instead it may indicate that variations involving the accre- 
tion disc (e.g., Wiita 2006) are important for them. We also 
concluded that there may be a weak negative correlation be- 
tween H/3 EW and I, but there is a significant one between 
the Mg II EW and I. In addition, we noted that there is a 
tendency for sources with detectable optical microvariabil- 
ity to have somewhat lower luminosities than those with no 
such detections (see Fig. 8 of Paper I) . 

All the above tentative conclusions from our Paper I 
were certainly interesting; however, as we stressed there, it 
is very important that they should be tested by examin- 
ing larger samples. The importance of a larger sample for 
such investigations is evident from our recent finding on the 
fraction of broad absorption line QSOs showing microvari- 
ability (Joshi et al. 2011). That fraction appeared to be 50 
per cent based on only 6 sources, but using a larger sample 
of 23 sources, Joshi et al. (2011) found it to be around ten 
percent, similar to RQQSOs (Gupta & Joshi 2005). 

Therefore, given the important implications of our 
above results, based on spectral analysis of 37 sources, for 
the origin of RQQSO microvariability, and hence to pro- 
cesses on the accretion disc itself, it becomes important to 
carry out spectral analyses of as many sources searched for 
microvariability as possible. The obvious source for such a 
sample would be the remainder of the objects from the to- 
tal of 117 sources in the compilation of Carini et al. (2007) 
for which SDSS DR-7 spectra are not available. This forms 
the main motivation of this paper, in which we more than 
double the sample size by taking new spectra and gather- 
ing others from the literature that we then analyze. Results 
from this larger sample can test the validity of results found 
using our previous modest sample in Paper L 

The paper is organized as follows. Section 2 describes 
the data sample and selection criteria. Section 3 describes 
the observation and data reduction while Section 4 gives de- 
tails of our spectral fitting procedure. In Section 5 we focus 
on BH mass measurements and in Section 6 we give esti- 
mates of Eddington ratios and of BH growth times. Section 
7 gives a discussion and conclusions. Throughout, we have 
used flat cosmology with Ho=70 kms~^Mpc^^, S7m=0.3 
and Qa=0.7. 



2 DATA SAMPLE AND SELECTION 
CRITERIA 

From the compilation of 117 radio-quiet AGNs in Carini 
et al. (2007), that were often extensively searched for mi- 
crovariability. Paper I analyzed 33 sources (in their total 
sample of 37 sources), for which SDSS DR7 spectra were 
available. Here we examined the rest of the the 84 sources, 
and using their redshifts found that among them 25 sources 
are such that neither their H/3 nor Mg ll lines fall in the 
optical spectral range as they have redshifts outside the fol- 



lowing respective ranges: z ^ 0.65; 0.43 z ^ 1.86. For 
the remaining 59 sources we found H/3 and Mg ll emission 
lines fall in the spectral range of 3800A-8300A for 50 and 9 
sources, respectively. 

We have searched for the optical spectra of these 59 
sources in SDSS DR8 (Aihara et al. 2011) and as well as 
in the NED data base. We could find SDSS DR8 spectra 
for 8 sources and NED spectra with desired quality (good 
S/N and spectral coverage that cover Mg ll or H/? line) for 
5 sources as listed in Table [1] 

We planned the observations of the rest of the 
46 sources, using the lUCAA Faint Object Spectro- 
graph (IFOSG) mounted on the 2-m telescope at lU- 
CAA Girawali Observatory (IGO), near Pune, India. 
With IGO we could obtain the desired spectra for 33 
sources. Among the remaining 13 sources: 6 were not 
visible from IGO, having very negative declinations; 4 
(J001555. 1+023024, J012017.2-H213346, J220311. 5-180143 
and J194240. 6— 101925) could not be observed either due 
to non-visibility or lack of enough observing time; 1, 
namely J103206. 2-1-324015, was classified as a white dwarf 
in NED; and the remaining 2, J082740.2+094208 and 
J154559. 1+270630, though observed, were dropped from our 
analysis as they showed very poor H/3 lines in their spectra. 

So finally we are left with a sample of 46 sources (33 
new observations, 8 from SDSS DR8 and 5 from NED); 
among them, 42 spectra cover H/3 lines while 4 spectra cover 
Mg II doublet lines. The source information is listed in Ta- 
ble [1] The first four columns give source name, RA(J2000), 
DEC(J2000) and emission redshift (zem), while the fifth col- 
umn provides the apparent magnitude in either B or V band. 
The sixth column provides the absolute B magnitudes and 
the seventh column indicate microvariability detection sta- 
tus; the eighth column gives a source classification from 
among classes such as QSO, BALQSO (Broad Absorption 
Line QSO), or Seyfert (Sy) galaxy type. In the last column 
we have given the spectral resource information. 



3 OBSERVATIONS AND DATA REDUCTION 

The new spectra were obtained for 33 sources. The observa- 
tions were carried out using the IFOSC mounted on the 2 
meter IGO telescope. We took long-slit spectra covering the 
wavelength range 3800 - 6840 A using the GR7 grism and 
5800 - 8300 A using the GR8 grism of the IFOSC to cover 
the Mg II and H/? lines, respectively, as noted in the third 
column of the observing log in Table [2] A slit width of ei- 
ther l."0 or l."5 was used, as listed in the fourth column of 
Table [21 Typical seeing during our observations was around 
l."2 to l."3. The raw CCD frames were cleaned using stan- 
dard irafQ procedures. We used halogen flats for flat field- 
ing the frames. Since at A > 7000A simple flat fielding does 
not remove the fringes, therefore in case of grism GR8 the 
QSO was moved along the slit for two different exposures at 
the same position angle. We subsequently removed fringing 
by subtracting one frame from the other taken on the same 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, wliicli are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. New sample of radio-quiet QSOs and Seyfert galaxies from the compilation of Carini et al. (2007) 
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" Object name in J2000 

Redshift as determined from the peak of O IIIA5007 emission line. 
* Apparent magnitudes (mg or my) and absolute magnitudes (Mg) are taken from Carini et al. (2007). 
<^ NED spectral references: ^Moustakas & Kennicutt (2006); ^'^'-fRoroson & Green (1992) ; *Kim et al. 
(1995). 
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night. The same procedure was applied to standard stars as 
well. We then extracted the one dimensional spectrum from 
individual frames using the IRAF task "doslit" . Wavelength 
calibration of the spectra was performed using Helium-Neon 
lamps. The spectrophotometric flux calibration was done us- 
ing standard stars and assuming a mean extinction for the 
IGO site. In cases when multiple exposures were needed we 
coadded the flux with weightage, where ai is the er- 

ror on the individual pixel. The error spectrum was com- 
puted taking into account proper error propagation during 
the combining process. The typical median signal-to-noise 
(SNR) over the spectral fitting range of our sample vary 
from 20-150 per pixel, as listed in last column of Table [2] 
Our spectral range also covers atmospheric absorption line 
regions but in all cases our fitting regions fall either bule- 
ward or redward of them, so we need not correct for any 
blending due to them in our analyses. 



4 LINE ANALYSIS THROUGH 

SIMULTANEOUS SPECTRAL FITTING 

The H/? and Mg ll emission lines were carefully modelled by 
using the fitting procedure we used in Paper I. The spectra 
were first corrected for Galactic extinction using the extinc- 
tion map of Schlegel et al. (1998) and the reddening curve 
of Fitzpatrick (1999). Then they were transformed into the 
rest frame using the redshift as determined from the peak 
of O IIIA5007 emission line. Limited by the complications 
needed to fit the continuum and the Fe ll emission, viz. (i) 
there are essentially no emission-line free regions where the 
continuum can be determined (Vanden Berk et al. 2001); (ii) 
the prominence of Fe ll features and their blending with the 
H/3 and Mg ll lines; (iii) the H/3 line is highly blended with 
the O lllAA4959,5007 lines, we have opted to carry out simul- 
taneous fitflof continuum, Fe ll emission, H/3 and Mg ll and 
all other metal emission lines present in the spectra. For this 
purpose we adopted the procedure as described in detail in 
Paper I, which in brief is as follows. 

4.1 H/3 region fit 

(i) We fit the spectrum comprising the H/3 line in 
the rest wavelength range between 4435 and 5535 A. The 
continuum in this region is modeled by a single power law, 
i.e., aiA~". 

(ii) The complex profile of the H/3 line is fitted with 
multiple (one to four) Gaussian with an initial guess of 
two narrow, one broad and one very broad components. To 
reduce the arbitrariness of the component fits and also to 
make the decomposition more physical, we have constrained 
the redshift and width of the two narrow components of 
H/3 to be the same as those of the O lllAA4959,5007 lines. 
The line profile of O lllA5007 (and hence of O lllA4959 ) 

^ To carry out the simultaneous fit we have used the 
MPFIT package for nonfinear fitting, written in In- 
teractive Data Language routines. MPFIT is kindly 
provided by Craig B. Markwardt and is available at 
|http://cow.physics.wisc.edu/'craigm/idl/| 



Table 2. IGO Observation log. 
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are modeled as double Gaussians, with one stronger narrow 
component with width less than 2000 km s~^ and one 
weaker broader component with width less than 4000 
km s~^ . However, if for some source spectra the second 
component is not statistically required, the procedure we 
use automatically drops it during the fit. So, we have 
only four free parameters, two for redshift and two for 
width, of the O lllAA4959,5007 and narrow H/J components. 

(iii) The optical Fe ll emission is modeled as C(A) = 
Cf,Cf,(A) -I- c„C„(A), where Cf,(A) represents the broad 
Fe II lines and C„(A) the narrow Fe ll lines, with the rel- 
ative intensities fixed at those of IZWl, as given in Ta- 
bles A.l and A. 2 of Veron-Cetty, Joly & Veron (2004). The 
redshift of the broad Fe ll lines and the H/3 component 
is fitted as a free parameter, while their widths are kept 
the same, with a constraint that they should be larger than 
1000 km s~^ ; this width is then used for estimating the BH 
masses. Similarly, for the narrow Fe ll line the redshift is 
fitted as a free parameter but the width is kept the same 
as that of the stronger narrow O lllA5007 component. The 
fourth (very broad) H/3 component, if required for the fit, is 
subject only to the constraint that its width should be more 
than 1000 km s~^ . The emission lines other than Fe ll and 
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Figure 1. The best fit to the H/3 emission line of the IGO spectra of the QSO J141348. 3+440014. Upper panel: complete spectrum fit 
(thick dashed/red) and components of the fit: power law continuum (thin short-dashod/black line), broad Fe II (thin long-dashcd/bluc), 
narrow Fe II (dotted/magenta) lines. Lower panel: continuum, Fe II and metal line subtracted spectrum (solid/black) with the best fit 
total H/J profile (thick dashed/red) and H/3 components (dotted/blue) lines. Note that the entire fit is performed simultaneously (not 
first continuum subtraction then H/3 fit) but these aspects are shown separately for the sake of clarity. 



O lllAA4959,5007 (see Table 2 in Vanden Berk et al. 2001), 
are modeled with single Gaussians. 

The final fit is achieved by simultaneously varying all 
the free parameters to minimize the value until the re- 
duced X? is ~ 1. Samples of our spectral fitting in the optical 
region are given for H/3 in Fig. [1] and Fig. [2] The values for 
FWHM and EW (both the broad component, EWb, and 
total, EWaii) for the H/? lines are given in Table O 



4.2 Mg II doublet region fit 

(i) We fit the spectrum comprising the Mg ll doublet 
region in the rest wavelength range between 2200 and 3200 
A. The continuum in this region is modeled by a single 
power law, i.e., aiA^". 

(ii) For fitting the Mg ll doublet we have used a Gaus- 
sian profile model (Salviander et al. 2007) with the ini- 
tial guess of two Gaussian components for each line of the 
Mg II doublet; however, if the second component is not 
statistically required the procedure automatically drops it 
during the fit. The redshift and width of each component 



(narrow/broad) of Mg llA2796 were tied to the respective 
components of the Mg llA2803 line. The peak intensity of 
Mg IIA2796 was constrained to be twice that of Mg llA2803 
as is predicted theoretically. 

In addition, we have constrained the width of narrow 
component to be smaller than 1000 km s~^ and the width 
of the broader Mg ll component to be same as width of UV 
Fe II emission line in the region (UV Fe ll ). We used an UV 
Fe II template generated by Tsuzuki et al. (2006), basically 
from the measurements of I ZW 1, which also employ calcula- 
tions with the CLOUDY photoionization code (Ferland et al. 
1998) . This template is scaled and convolved to the FWHM 
value equivalent to the broad components of Mg ll by tak- 
ing into account the FWHM of the I ZW 1 template. The 
best fit value of the broad component of Mg ll obtained in 
this way is finally used in our calculation of BH mass. To 
test for any overfitting caused by assuming two components, 
we also forced our procedure to fit only single components; 
however, in doing so for all our sources, a very good fit is 
never found for the wings of the lines nor for their central 
narrow cores. 

(iii) Emission lines other than Fe ll lines identified from 
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Figure 2. As in Fig.[T]IGO spectra for J153638.3+543333. 



the composite SDSS QSO spectrum (see Table 2 in Vanden 
Berk et al. 2001), are modeled with single Gaussians. 

The final fit is achieved by varying all the free param- 
eters simultaneously by minimizing the value, until the 
reduced Xr is ~ 1- Demonstrations of our spectral fitting in 
the UV region are given in Fig. [S] and values for FWMHs 
and EWs are given in Tabled 

4.3 Optical microvariability and spectral 
properties 

The results from Paper I, with a modest sample of 37 
sources, showed that the spectral properties for the sources 
with and without optical microvariability are quite similar. 
Here we search for better statistical results by adding our 
new 46 sources to the sample in Paper I. In these added 46 
sources, among the 42 sources for which we have spectral 
coverage of the H/3 line, 9 have shown optical microvariabil- 
ity while the other 33 have not been seen to show this prop- 
erty (Table [1]). Of the 4 new sources with spectral coverage 
of the Mg II doublet, not one showed optical microvariability. 
We now investigate the correlation between spectral proper- 
ties (i.e., FWHM and EW) and optical microvariability with 
our best— fitting values and with the larger combined sample 
of 83 sources. Among the 53 sources with only H/3 line cov- 
erage, optical microvariability is shown by 15 sources, while 



of the 22 sources with only Mg ll doublet spectral coverage, 
only 4 have been shown to exhibit microvariability prop- 
erties. Of the remaining 8 source spectra that cover both 
Mg II and H/3 lines, optical microvariability was shown by 
one source. 

From the multiple component fit of the H/3 line we have 
used the broad component fit (Sect. I4.1|l to perform the 
comparison. This is because the clouds responsible for this 
broad component are clearly in the sphere of influence of 
the massive BH, while the other components may not be. 
In Fig. |4] we show the histograms of FWHM and rest frame 
EW values based on our best fits for H/3 and Mg ll lines. The 
shaded and non-shaded regions correspond respectively to 
sources with and without confirmed optical microvariability. 
From these plots it appears that the distributions of sources 
with and without microvariability are on the whole quite 
similar. To quantify any differences in these distributions 
we have performed Kolmogorov-Smirnov (KS) tests on all 
the distributions shown in Fig. 2] For the null hypothesis 
that the samples are drawn from the similar distributions 
we found the probabilities for the two H/3(EW) distributions 
to be 0.90 and for the two Mg II(EW) distributions to be 
0.70. Similarly, we found the KS-tests null probability value 
for the two H/3(FWHM) distributions to be 0.29 and for the 
two Mg II(FWHM) distributions is 0.99. 

In addition in the "blazar component" scenario one 
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Figure 3. Mg II emission line fits for two RQQSOs, J125948.8+342323.5 and J125048. 3+395139. The upper plot in each panel shows the 
SDSS spectra and continuum as a dashed line with the best fit UV Fe II template as a solid (magenta) line. The lower plot in each panel 
show the continuum and Fe II subtracted spectrum (sohd) the best fit Mg II profile (thick dashed/red) and Mg llAA2796,2803 components 
by (dotted/blue) fines. Note that the entire fits are performed simufiancously (not first continuum subtraction, then Mg II fit) but are 
shown separately for the sake of clarity. 



would also expect some dilution of emission line strength 
of the Fe ll template for microvariable sources. To test this 
possibility we have shown in Fig. [5] the EW distribution of 
broad Fe ll emission template for variable and non-variable 
sources, which by eye appear indistinguishable. The KS-test 
null probability value for these two Fe ll distributions is 
0.72, and so does not give any hint of a blazar component. 
We have also tried to investigate the difference in fraction of 
broad components of H/3 and Mg ll to the respective total 
line intensity for the sources with and without microvari- 
ability. We found that the distribution of this fraction of 
the H/3 and Mg ll broad component is quite similar both for 
variable and nonvariable sources as is shown in Fig. [S] 

The above substantial null-hypothesis probabilities do 
not establish any relation of EW or FWHM with the optical 
microvariability properties and thus do not support a blazar 
component model of microvariability for RQQSOs, which 
would predict smaller EW values for variable sources, due to 
dilution of emission line strength by jet components (Czerny 
et al. 2008). Here, the results drawn from a sample more 
than twice as large are in agreement with our results in 
Paper I. 



5 BLACK HOLE MASS MEASUREMENTS, 
EDDINGTON RATIOS AND BLACK HOLE 
GROWTH TIMES 

We have estimated the black hole masses using the virial 
single epoch method (Dibai 1980), following the approach in 
Paper I. This method has been shown to be consistent with 
reverberation mapping masses (e.g., Bochkarev & Gaskell 
2009). Improved empirical relationships (e.g., Vestergaard & 
Peterson 2006) are used for estimating the central BH mass 
based on FWHMs of emission lines such as H/3 and Mg ll. 
We used Eq. 1 and Eq. 2 in Paper I, taken from Vestergaard 

6 Peterson (2006) and McLure & Dunlop (2004), to deter- 
mine respectively black hole masses based on H/3 lines and 
Mg 11 lines, viz; 



log Mbh(H/3) =log 



/ FWHM(H/3) \ ' 
i 1000 km s-i ) 



/ , r ! N \ 0.50 + 0.06 

+ (6.91 ±0.02) + log ''f'''^ ) 
\ lO-i^erg s-i / 



(1) 
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Table 3. Results of the H/3 line analysis 
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Table 4. Results of the Mg II line analysis 



QSOName z,^, loXT^-i FWHM(kms-i) Log(^)^ ^^si^t T± EWb(A) EW„n{A) 



J112731.9-304446 


0.6679 


5691.58 


1447.36 


8.53 


9.00 


2.32 


13.14 


13.14 + 0.63 


J125048. 3+395139 


1.0318 


1294.13 


3029.19 


8.78 


9.32 


0.25 


19.64 


23.46 + 0.31 


J125948.8+342323 


1.3762 


2577.21 


3297.04 


9.04 


9.54 


0.30 


14.12 


16.28 + 0.27 


J135458. 7+005211 


1.1253 


1451.56 


3007.11 


8.80 


9.34 


0.27 


22.30 


26.21 + 0.17 



Using the McLure & Dunlop (2004) scaling relation, i.e., Eq. [2] 
^ Using the fixed slope of the i — L relation from Dietrich et al. (2009), their Eq. (6). 
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Figure 4. Histograms of FWHM and rest frame equivalent width (EW) based on best fits of H/3 and Mg II lines. The red (dark) shaded 
regions correspond to sources with confirmed optical microvariability while the striped regions correspond to those sources for which 
optical microvariability has not been detected. 
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Figure 5. Histograms of EW of broad Fe II template . The red 
(dark) shaded regions correspond to sources with confirmed opti- 
cal microvariability while the striped regions correspond to those 
sources for which optical microvariability has not been detected. 
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Figure 6. The left panel shows the histograms of the fraction of 
the broad H/3 component in the H/3 emission lines and the right 
panel shows the broad Mg II component fraction in Mg II emission 
lines. 
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Figure 7. The upper left and middle panels show the histograms of Eddington ratio L^j^i/ L^dd a-nd logC^^bh/A/©) based on the H/3 line, 
while the corresponding bottom panels show the results based on the Mg II line. The shaded regions correspond to sources with confirmed 
optical microvariability while the non-shaded regions are for those sources for which optical microvariability is not detected. The right 
upper and bottom panels respectively show plots of log^Mf^i^ / Mq) versus I/joi. Triangles and circles respectively show sources with and 
without confirmed optical microvariability. 



log AfBH(Mg II) = (0.62 ± 0.14) log 



f AZ,a(3000A) 
I lO-i^erg s-i 



where La (5100 A) & La (3000 A) are the monochromatic lu- 
minosity at 5100 A and 3000 A , respectively, which we have 
computed from the best fit power-law continuum, a^X'" , in 
our simultaneous fit of the whole spectral region. 

We have also estimated the Eddington ratio I = 
Lboi/Ledd, where Lboi is taken as 5.9 x ALa(3000A) and 
9.8 X ALa(5100A) for Mg ll and , respectively (McLure 
& Dunlop 2004), and Ledd = 1.45 x 10^*(M6h/Afo) erg s'S 
assuming a mixture of hydrogen and helium so the mean 
molecular weight is fj, = 1.15. The combined sample results 
are given in Fig. [7] which by eye indicates that distribu- 
tions of sources with and without optical microvariability 
appear similar with respect to both BH mass and £. Quanti- 
tatively, this is also supported by KS tests, which show that 
the probability of the null hypothesis (Pnuii) for sources with 
and without microvariation, is as high as 0.69 for BH mass 
distributions and 0.11 for Eddington ratio {£) distributions 
based on H/3 lines fit. Similarly using Mg ll lines, the null 
hypothesis (Pnuii) for sources with and without microvaria- 
tion, is as high as 0.98 for BH mass distributions and 0.94 
for Eddington ratio {£) distributions. 

To test the reasonableness of estimated Eddington ra- 



tios, we also computed black hole growth times to compare 
them with the age of the Universe (at the time the QSO is 
observed), by using the following equation (Dietrich et al. 
2009), 



Mtnit^ts) = M^nto) exp ( I ^^—^ 

e tedd 



(3) 



where r = tots — to is the time elapsed since the initial 
time, to, to the observed time, tabs', M^f^"^"^ is the seed BH 
mass; e is the efficiency of converting mass to energy in the 
accretion flow, and t^dd is the Eddington time scale, with 
tedd = aTc/4nGmp = 3.92 x 10** yr (Rees 1984). We used 
Eq. |3] to derive the times, r, necessary to accumulate the 
BH masses listed in Tables |3] and U for seed black holes 
with masses of M^^""^ = IOMq, 10^ M©, and 10^ Mq, re- 
spectively. Two cases are considered: (i) BHs are accreting 
at the Eddington- limit, i.e., £ = Lboi/i^edd = 1-0 and the effi- 
ciency of converting mass into energy is e = 0.1; (ii) BHs are 
accreting with our observed Eddington ratios and e = 0.1. 
These results are summarized in Tables [5] and [S] showing 
that the key criterion that black hole growth time should be 
smaller than the age of Universe is fulfilled for all sources. 

In our Paper I we found a weak negative correlation 
between H/3 EW and the Eddington ratio, £, and a signif- 
icant one between the Mg ll EW and £. In Fig. |8] we show 
the observed variations of FWHM and EW with £ based on 
both H/3 and Mg ll lines for our combined sample. There ap- 
pear to be linear relations of both FWHM and EW with £ in 
these log-log plots. For the FWHM plots this is unsurprising 
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Values of I estimated using H/3 lines (Table [Sjl. 



since the BH masses are proportional to L^dd- To quantify 
any such linear relations for the EW plots we perform linear 
regressions, treating I as the independent variable, and find 

logEW(H/J) = (1.33 ± 0.12) + (-0.05 ± 0.13)log(^), 
logEW(Mg II) = (1.07 ± 0.04) + (-0.29 ± O.lO)log(^). (4) 

Here the errors on the fit parameters are purely statisti- 
cal. We have calculated the Spearman rank correlations of 
logEW with log^, and found the correlation coefficient for 
Ts{H/3) = —0.05, with null probability Pnuii ~ 0.72, and so 
no significant correlation is present. Whereas rs(Mg II) = 
—0.45, with pnuii = 0.013 and so this negative correlation is 



significant. These results are in agreement with those from 
our Paper I that were based on a modest sample size. 



6 DISCUSSION AND CONCLUSIONS 

Modern surveys such as SDSS have allowed investigators to 
carry out spectral analyses of large numbers of quasars to es- 
timate central BH masses using virial approaches (e.g., Shen 
et al. 2008; Fine et al. 2008; Vestergaard & Osmer 2009) 
and to understand their demography (Dong et al. 2009b). 
These studies have resulted in very important insights on 
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Table 6. BH growth time estimates (in units of Gyr) from the Mg II hne analysis 
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Figure 8. Observed variation of FWHM (upper) and EW (lower) with Eddington ratio {£ = ^ibol/^edd) based on both H/3 (left) and 
Mg II (right) lines. Triangles and circles respectively show sources with and without confirmed optical microvariability. The dashed lines 
show the linear regression fits treating £ as the independent variable. The correlations of both FWHM and EW with £, suggest that the 
Eddington ratio is one of the fundamental parameters responsible for some AGN properties. 



the important physical parameters of AGN central engines 
and their environments. 

Dong et al. (2009b) found that the variation of the 
emission-line strength in AGNs are regulated by £, presum- 
ably because it governs the global distributions of the prop- 
erties such as column density of the clouds gravitationally 
bound in line emitting region. Shen et al. (2008) have found 
that the line widths of H/3 and Mg 11 follow lognormal dis- 
tributions with very weak dependencies on redshift and lu- 
minosity. Fine et al. (2008) used Mg 11 lines to estimate BH 



masses, and found that the scatter in measured BH masses 
is luminosity dependent, showing less scatter for more lumi- 
nous objects. 

The sample we have considered more than doubles our 
sample in Paper I, but still is much smaller compared to 
those in the above papers, that consist of between 1100 and 
almost 57,700 quasars. However, each member in our sample 
has been carefully selected to be among the special group 
of RQQSOs and Seyfert galaxies already examined for opti- 
cal microvariability (e.g., Carini et al. 2007). This criterion 
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Figure 9. As in Fig. |8] with Lbol as the independent variable. This panel indicates that the sources with and without optical microvari- 
ability are similar in luminosity. 



demands that modest aperture (usually 1-2 m) telescopes 
can make precise photometric measurements in just a few 
minutes, and so limits the members to the rare QSOs with 
bright apparent magnitudes (usually my < 17.5). In addi- 
tion, we have taken special care in the fitting of the line 
profiles as discussed in Sections 4.1 and 4.2. 

It has been found that for BL Lacertae objects the emis- 
sion line detection and strength varies with overall contin- 
uum flux (e.g., Nilsson et al. 2008). For instance, when BL 
Lacs are in optically faint states weak emission lines have 
sometimes been detected that are usually not seen during 
their high states, presumably due to their being swamped 
by the Doppler boosted continuum arising in a strong jet 
component (Nilsson et al. 2009). 

To check whether the continuum states of our ob- 
jects, high (bright) or low (faint), have any such effect on 
their spectral properties, we have searched in the SDSS 
DR8 archive for multi-epoch photometric and spectroscopic 
observations for the 44 SDSS sources in our sample. We 
found multi-epoch photometric fluxes were available for 14 
sources, with time gaps ranging from a few days to years. We 
then searched for spectroscopic observations of these sources 
and found multi-epoch spectra are available only for three 
sources, viz. J025937.46+003736.3; J084030.0-H465113; 
J093502. 6-1-433 111. Among them, J093502. 6-^4331 11 has 
two spectra with a roughly one year time gap (observations 



on 2002 Feb 20 and 2003 Jan 31) but only a single photo- 
metric data point at all close to either of those dates (2001 
Dec 20). For J084030. 0-1-465113 two spectra were available 
about a month apart (2001 Jan 13 and 2001 Feb 19) but 
these had their two closest photometric observations made 
on 2001 Jan 24 and 2001 Jan 26, respectively. As these two 
photometric observations have only a two day time gap and 
have a g-mag difference of only 0.01 this object also does 
not serve to test for such a correlation. We are left with just 
one source, J025937.46-I-003736.3, that has both photomet- 
ric and spectroscopic multi-epoch data that are well paired. 
For this source there were two spectral observations on 2000 
Nov 25 and 2001 Sep 21 and corresponding photometric ob- 
servations on 2000 Nov 27 and 2001 Sep 21 with g— mags of 
16.41 and 16.47, respectively. As this g— mag difference of 
0.06 is much larger than the typical RQQSOs INOV mag- 
nitude variation of about 0.01 mag over a night, this source 
is the only one that might test whether brighter or fainter 
states have any significant effect on spectral properties. Even 
though this difference in brightness is modest, we carried 
out our spectral fits to each of those two spectra to esti- 
mate the FWHMs and EWs of the H/3 line, using the same 
method as we used for other sources (Sect. 4.1). The best fit 
EW(H/3 ) of broad component for first epoch was found to 
be 30.47 ±0.63A and for second epoch, 29.23 ±0.45A. These 
values are statistically indistinguishable, so there is no ev- 
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idence for any effect due to variation in source brightness. 
However, in order to say anything firm about this possibihty 
for RQQSOs, nearly simultaneous spectroscopic and photo- 
metric measurements would need to be made on at least two 
occasions for a decent sized sample. 

In this paper we continued our work begun in Paper 
I with 37 RQQSOs by analyzing spectra for an additional 
46 sources. Most of these spectra were obtained by us at 
the IGO. We conclude that there is a significant negative 
correlation between Mg ll EW and t, (Fig. El Section 5), 
as also found in Paper I (there Fig. 7) and by Dong et al. 
(2009b). However, we have not found any significant corre- 
lation between the equivalent widths of the H/3 lines and the 
Eddington ratio, I. We can also see from Fig. [8] that there is 
a decline in FWHM with I; this is not surprising since the 
BH masses are proportional to Ledd- In Paper I (Fig. 8), we 
noticed a interesting tendency for sources with detectable 
optical microvariability to have somewhat lower luminosity 
than those with no such detections, which required investi- 
gation with a larger sample. One might expect such a trend, 
as lower mass BHs would have correspondingly shorter phys- 
ical timescales. However, as can be seen from Fig.[9]there is 
no such trend seen with our new larger sample of microvari- 
able sources (16 as compared to 7 in Paper I). 

We also find that the BH masses estimated from the 
FWHMs of both the H/3 and Mg ll lines are reasonable, in 
that growth to their estimated masses from even small seed 
BHs are easily possible within the age of the Universe at 
their observed redshift if the mean I values are close to unity 
(Tables 5 and 6). This remains true for the great majority of 
RQQSOs even if the value of I we compute from the current 
continuum flux was constant until the time we observe them; 
however, this assumption does not work for 8 out of the 42 
QSOs with H/3 lines, while we consider Mbh{seed) — IOA'Iq, 
for 4 while we consider Mthiseed) = 10^ M0 and for none 
while we consider Mi,h{seed) = 10^ Mq, suggesting that in 
a few cases Mbh{seed) could need to be as large as IQp Mq 
or the accretion rate was substantially higher in the past. 
With Mg 11 lines profiles no QSO was found problematic 
with respect to this assumption (Table 6). 

As Fig. 13] shows, histograms of FWHM and rest frame 
EW values for sources with and without confirmed opti- 
cal microvariability are on the whole quite similar. Un- 
der the null hypothesis that the samples are drawn from 
the similar distribution using KS-tests we find a probabil- 
ity of 0.90 for the H/3(EW) distributions and 0.70 for the 
Mg II(EW) distributions. Similarly this null probability for 
the FWHM distributions is 0.29 for H/3(FWHM) and 0.99 
for Mg II(FWHM) distributions. We conclude that EW or 
FWHM distributions (for both the H/3 or Mg 11 ) are prob- 
ably independent of the presence or absence of detected mi- 
crovariability properties in those QSOs and Seyferts. As dis- 
cussed in the Introduction and in more detail in Paper I, if 
much of the optical emission in RLQSOs comes from a jet, 
then we would expect the EWs of the RLQSOs to be sig- 
nificantly lower than those of the RQQSOs and that the 
EWs of microvariable sources would be less than those of 
non- variable sources. Our results are in agreement to the 
conclusion in Paper I and thus do not support the hypothe- 
sis (e.g., Gopal-Krishna et al. 2003; Czerny et al. 2008) that 
RQQSOs possess jets that are producing rapid variations. 



Instead it may indicate that variations involving the accre- 
tion disc (e.g., Wiita 2006) play an important role here. 

Further improvements to our results could be obtained 
through extensive searches for INOV in a larger sample of 
RQQSOs to reduce the statistical uncertainties. In any such 
studies it would be most useful to take the spectra just be- 
fore or after the photometric monitoring run, so as to ensure 
the simultaneity of light curve and spectral properties. This 
would rule out any possibility of change in either of these 
properties due to temporal gaps between the spectral and 
photometric epochs; while we very much doubt that this is 
an important effect, it could have an impact on our results. 
Such larger samples could be optimally designed if they were 
made as homogeneous as possible on the basis of apparent 
magnitudes, redshifts and absolute magnitudes. 
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